function y = culD(x,H,m)
    a = 2.518892
    b = 0.2980319
    c = 3.059665
    y = b.*H.^2+c.*exp(m)+a
end